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NASA TECHNICAL MEMORANDUM X- 64651 


GRAVITY GRADIENT TORQUE PROFILES OVER AN ORBIT FOR 
ARBITRARY MODULAR SPACE STATION CONFIGURATIONS IN 
THE Y-POP AND INERTIAL HOLD MODES 

INTRODUCTION 


The modular space station is currently receiving a great deal of study 
by NASA^» ^ and its cdntrol is obviously a topic of considerable importance. 

A major contributor to the control problem is the gravity gradient torque, 
which is considerable for a large orbiting body and which, for that matter, ’can 
cause instability of even very small satellites if they are improperly oriented 
[ 1 , 2 ] . 

Prevention of the continual acquisition of angular momentum caused' by 
such torques clearly requires some advance knowledge of the magnitudes of the 
torques to be controlled. The objective of this report is to examine the gravity 
gradient torque distribution acting on an arbitrary configuration of the modular 
space station over any orbit. This is done for both the Y-POP configuration 
(in which the space station moves with its longitudinal X-body axis always 
directed along the velocity tangent and the Z-body axis along local vertical) and 
for the station frozen into and remaining in any inertial hold configuration during 
the above-mentioned Y-POP orbit. 


The procedure used and described here allows for solar panel gimbaling 
and for varying the values of the spatial and orbital parameters F^, 

and orbiting angle 9 . Any Initial Space Station (ISS) configuration may be 
built up, and with modifications the procedure can be made 'to 'accommodate the 
Orowth Space Station (GSS), as well as other orbital modes. 


V 

1. Space Station First Performance Review. MSFC-DPD-235/DR NO. MA-02, 
MDC G2279, McDonnell-Douglas Astronautics Co., April 1971. 

2. Space Station Mass Property Status Report. MSFC-DPD-235/DR NO. 

SE-07, MDC G2313 McDonnell-Douglas Astronautics Co. , June 1971. 



It is shown by examples that solar panel gimbaling can have a significant 
effect on gravity gradient torque distribution and hence on momentum require- 
ments. The examples also show wide variations in bias torque for various 
stations examined. 

OUTLINE OF THE GENERAL PROCEDURE 


Before the details of this study are discussed, the general procedure 
that was followed is outlined below for the convenience of the reader: 

1. Give or vary the angles F , 0 , \ , A, , and 0 of Figure 1. 

y z y z 

2. Compute the quantities (Ui, U 2 , U 3 ) from equations (21) through 

(23). 

3. Compute 17 ^ and (Fig. 1) from equations (15) and (16). 

X y 

4. Compute the gimbal angles 61 and 63 (Fig. 2) which the solar 
panels must turn through to face the sun [equations (30) and (33)] . 

5. Transform the solar panel inertia tensor to vehicle coordinates 
(Fig. 2). 

6 . Given the configuration of the station, build it up, determine its 
mass center C , and compute its composite inertia tensor with respect to C . 
Perform this in two ways as a check. 

7. If the vehicle is in inertial hold, transform the inertial properties 
to the orbiting coordinate system T (Fig. 3). 

8 . Compute the principal moments of inertia and their directions, 
using a standard library routine. 

9. Compute the gravity gradient torques acting on the vehicle. 

10. Increment 0 (orbiting angle) and (to account for orbital 

regression) and repeat steps 7 and 9 for inertial hold and all but step 7 for 
Y-POP, Continue for a complete orbit, and repeat the entire procedure for 
various configurations and parameters of interest. (In the Y-POP mode, only 
the solar panels* changing contribution to the inertia tensor is recomputed in 
step 6 . ) 
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VARYING INERTIAL PROPERTIES OF THE SOLAR PANELS 


This section describes the gimbaling of the solar panels necessary in 
order that they always face the sun. Once these angles are computed, they are 
used to transform the solar panel inertia tensor to the vehicle axes. 

The two gimbal angles for the panels are shown in Figure 2 . The first 
angle is a rotation about X of amount 6i . Note that (X , Y , Z ) are 

oJr olr oP oP^ 

axes parallel to the vehicle axes (X^, Y^, Z^) , passing through the center of 

mass of the panels. This first rotation brings the panel axes (frame SP) into 
the primed axes (frame B' ) shown at the upper right in Figure 2 . The second 
rotation is about the Y' axis, through the angle 62 > bringing the solar panel 
symmetry axes into the (X" , Y" , Z" ) positions (frame B" ), shown at the 
lower right of Figure 2 . 

Since the inertia tensor with respect to reference frame B" is known, 
we can thus transform these properties back through B' to SP and finally to 
the V-frame, where all the inertias are being collected. 

The condition which must be fulfilled here is that unit vector id’ (in the 
Z” direction) is always pointing toward the sun; i. e. , it is equal ^ > which 

is defined to be a unit vector in the Z direction. This then gives two addi- 

XV 

tional equations stating that the components of ic" in the X and Y direc- 

R R 

tions be zero at all times. From these equations we will arrive at a solution for 
the gimbal angles 61 and 62 > as will be shown shortly. 
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We are thus ultimately seeking a transformation matrix [ ]VLg„ which 

will express the vinit vectors of B" in terms of those of the solar reference 
frame (X , Y„, Z„), which is labeled R . In steps, we have 

Iv R R' 




( 1 ) 


in which subscript O refers to the orbital plane frame (Fig. 1 ), having coord- 
inates (X , Y , Z ). 

- ' o o o' 

Let us first compute matrix . Clearly, we have 


and inspection shows that 


( 2 ) 


and 


^^B" B'^ 




f“vo' = 


'cos 62 

0 

-sin 62' 

0 

1 

0 

sin 62 

0 

cos 62 

1 

0 

0 

0 

cos 61 

sin 61 

0 

-sin 61 

cos 61 


cos 0 0 -sin 0 

0 10 
sin 0 0 cos 0 


( 3 ) 


( 4 ) 


( 5 ) 


in which 0 is the orbiting angle shown in Figure 1 , which locates the space 
station in the orbital plane with respect to the ascending node. 


Substituting equations ( 3 ) through ( 5 ) into equation ( 2 ) yields 
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[M 


B ”0 


C;: C -S. C. 

62 0 O2 61 0 

61 62 


s. s 


S„ c 

61 0 

61 

61 9 

S, C +C,, c„ s„ 

-s„ 

-S, S„ +C„ C, C 

62 d 62 0 

61 62 

62 9 61 62 0 


in which C- means cos 62 » S means sin 0 , etc. 

O 2 0 

The other matrix factor of equation (1) is [M ] . Rotating in the orbital 

OIv 

plane about through the angle , and then rotating through about the 
rotated position of the axis which began coincident with brings the Z-orbital 
axis into alignment with the solar reference axis Z . Therefore, 










c s 

n T] 

ly >: 



( 7 ) 


Substituting equations (6) and (?) into equation (l) yields the desired 
result for [ . We need only the third row here which, written in equa- 

tion form, is 


k" = + Ajjjj 



( 8 ) 


in which the A. are functions of 61, 
, it follows that Aj = A2 = 0 and 
respectively become as follows: 


62. T)y> Slid 0 • If k" is to equal 
A3 = 1 . Written out, these three equations 
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Sr c + C. C- C + s. S - S C- C. C = 0 , '(9) 

62 9 TJy 61 62 0 T)y 62 e TJy 7?y 61 62 9 

s. o s s + C- c. s s - S- c. c - S. s n s 

62 9 T}^ T]y 61 62 9 T]^ 7 ]y 61 62 T]^ 62 9 T]y TJ^ 


+ c C S =0 

61 62 9 Vy 


( 10 ) 

t . 


and 


S. CS C +C.C. SS C +S. C_S -S. SC C 
62 9 7 ]y 77^ 61 62 9 T)y r)^ 61 62 7 )^ 62 9 7 J^ TJy 


+ C- C- C„C C =1 
61 62 9 T)^ TJy 


( 11 ) 


Even though equations ( 9 ) through ( 11 ) are unwieldy, they can be solved 

for 61 and 62 once we know rj and rj in terms of the parameters r , » 

^ y y 

\ and X . The value of ff will be incremented over an orbit but will remain 
Y 2 

constant through each incremental step, so it will be known whenever equations 
( 9 ) through ( 11 ) need to be solved. . ; 

Now consider Figure 4 . Letting (Ui, U2, U3) be the projections of k^^ 

* L 

onto the axes of the orbital reference frame, it follows that 


Ui = Sint? cos 77 
'y 'x 


U2 “ -sin 7]^ 


( 12 ) 

( 13 ) 


and 


U3 - cos 7]y cos T7^ 


( 14 ) 


Inverting these yields 


77 = tan 
'x 


-1 


-U, 


(u 7 +U 3 ^ 


( 15 ) 


7 



0 




m.iifr/'oe 


Figure 4. Relation between 


^ -foltis seen from equations (15) and (16) that it remains only to solve 

for Ui, Uo, and Uq in terms of r , \ , and X to be able to solve 

^ ^ • y ■ z ■ y’ z 

equations (9) through (11). We do this via the transformation 


= '“rs"“se‘'“eo’ • 

Si) 

in which S refers to the ecliptic and E to the equatorial plane, as depicted 
in Figure 1. 

) The matrices of e.^uation (17) which must be multiplied are 
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C,C^ 

Ip r 
^ y 

CO 

-c,s 

-S c 
* 


Vr 

y 

0 



( 18 ) 


and 




C S 

^ <f> 




( 19 ) 


"*E0' 


c c 
\ \ 
y z 


-s c 

\ X 

y z 


-c s 

X X 

y z 


s s 

X X 

y z 


(20) 


fo equation (18), ip is the angle about Z that must be rotated to bring 

R 

the X-axis into the orbital plane where X„ ■ is defined. This angle is of no con- 
sequence, however, as we need only the third row of [ . After multiplying, 

then, the required components of become 




ti ivii , , .... 

= Ui = 'C ■ ^cii su a, • (21) 

^ F0X rA'-A ’ Jt ii vr-iaf* M ly 

y zyz yzz yy z 


(X^ - component) 




= Uo=-S CCS +s sc +c ss 
^2 r 0 ^ r X X 

y zyz yzz yyz 


( 22 ) 


(Y^ - component) 


and 




= Ug = s„ c, s, + c„ c. 


r ^ X r X 
y z y y y 


(23) 


( - component) 


Thus we may determine; (1) the U. from equations (21) through (23), 

(2) 7) and rj from equations (15) and (16), and (3) 6j and Sg from 
X y 

equations (9) through (11). This last step goes as follows: Equations (9) 
through (11) may be represented by (with Sg= sin 6g . etc. ); 


SgKj + CjCgKg = 0 , 

S2K3 + C1C2K4 + SjCgKg = 0 

SgKg + CjCgKg + SjCgKg = 1 


(24) 

(25) 

(26) 


in which the are constants for any position being examined, listed in 
Table 1 for completeness, hi Table 1, A is 0-r] , which is the lead angle 

of the space station past the sunline’ s projection into the orbital plane. 

Equation (24) yields 

Sg = -CiCgKg , (27) 

and substituting equation (27) into equation (25) yields 

CiCgKio + SjCgKs =0 (28) 
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If cos 62 = 0 > it can be seen from equations (24) through (26) that 
gimbal angle 61 drops out of the equations. We are then left with the following 
resulting trivial cases (after evaluating various possibilities); 


T?x = Tjy = 0- I , and 62 = - I (31) 

or 

STT 7T 

’iy = — ^2 = 2 • (32) 


(30), 


Returning to the nontrivial case, equation (24) yields, with equation 


62 = tan ^ (- tan A cos 61 ) ± tt 


(33) 


There are thus four pairs (61, 62) given by equations (30) and (33), 
which satisfy equations (24) and (25). However, only two of these pairs also 
satisfy equation (26), the other two answers representing solar panels pointing 
away from the svm. 

The results of an investigation of the proper quadrants which should be 
used to avoid more than 180-deg traverse in either direction are shown in 
Table 2. 


TABLE 2. GIMBAL ANGLE QUADRANTS 


A 

Quadrant 

For n between 
'x 

0 and 7 t / 2 

For n between 
'x 

0 and -7 t / 2 

(8„ 6j) 

Quadrants 

(ai, 62) 
Quadrants 

1 

(1. 4) 

(4, 4) 

2 

(2, 4) 

(3. 4) 

3 

(2, 1) 

(3. 1) 

4- 

(1. 1) 

(4, 1) 
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It is now a simple matter to transform the inertial properties from the 
B" frame to the SP frame by the standard tensor transformation 


I.. 


SP 



3 


r 

i=i 



(34) 


in which a is the direction cosine of the angle between the X " axis and 
SP pq P 

the axis. The computer program performs the transformation given in 

equation (34). 


INERTIA TENSOR FOR THE SPACE STATION 


. The computation of the gravity gradient torque depends upon a knowledge 
of the inertial properties of the body referred to its mass center. Therefore, 
this section briefly describes how the inertia tensor is computed for the space 
station. 


The last section described how inertial properties of the solar panel 
could be obtained with respect to axes through its center of mass parallel to 
the vehicle axes. This can be done even more easily for those modviles whose 
inertial properties are known with respect to axes parallel to those of the 
vehicle; these are the modules attached to ports numbered 1 or 4 in Figure 5. 



V 



DOCIC 
PORT 4 



Figure 5. Dock port labeling convention. 
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For the modules attached to docking ports 2 or 3, one lateral axis is 
considered pointed parallel to . Then, Mohr* s circle is utilized to trans- 
form the inertial properties, known with respect to the other two axes, to axes 
through the mass center of the module parallel to and . This is done 

within the computer program; no transformations need be performed externally. 

At this stage, we then have the inertial properties of panels and all 
modules with respect to their individual centers of mass. One by one, these 
are then transferred to the origin of the vehicle coordinate frame V by 
standard transfer theorems, such as 


' "x X ’ <“> 

c c 

^ Sc Z 
c c 

I 

in which (X, Y, Z) are axes parallel to the mass center axes (X^, Y^, Z^); 

m is the body* s mass; d is the distance between X and X^ , and a and c 

are the X and Z coordinates of the origin of (X, Y, Z) with respect to 
c c 

(X^, Y^, Z^). In this application, (X, Y, Z) = (X^, Y^, Z^). 

At the same time that the inertias are being transferred, sums are being 
kept of the systems* mass and of the mass moments with respect to (X^, Y^, Z^) . 

When all modules have been covered in this way, the mass center of the 
station is then determined simply by dividing the mass moments by the total mass. 
Then a transfer of the composite inertial properties from the V origin to the 
mass center of the station yields the inertias needed for the gravity gradient 
computations. 

Before we proceed to these computations, however, the inertia tensor 
is checked by going back and transferring each panel and module* s inertial 
properties directly to the now-known mass center. This will furnish a check 
not only on the inertia tensor but also on the mass center position calculations. 
This follows from the fact that if the mass center is incorrectly computed, the 
original inertia transfer from the V origin to this erroneous mass center will 
be incorrect because the formula used therein requires correct mass center 
coordinates. 
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If the station is orbiting in the Y-POP mode, the solar panel inertial 
properties are transformed and transferred in the manner just described at.^ 
each station of the orbit. However, the modules' inertial properties are com- 
puted only once since they are fixed with respect to the vehicle. 

When an inertial hold mode is called for, the inertial properties of -the 
entire station are then frozen with respect to the vehicle axes. However, these 
combined properties must be transformed to the T -frame, or orbiting frame, 
of Figure 3, because it is with respect to T that gravity gradient torque com- 
putations become rather simple, as described in a later section. Hence, the 
computer is taught to perform this inertia transformation, which depends upon 
the orbiting angle 0 . 

PRINCIPAL MOMENTS OF INERTIA AND CORRESPONDING 

PRINCIPAL DIRECTIONS 

This part of the procedure is concerned with determining principal 
moments of inertia which, along with their corresponding principal directions, 
are recjuired for determining gravity gradient torques. 

• ' Beginning with an inertia tensor now computed with respect to axes 
through the space station’ s mass center, the computer program then calls on 
the UNIVAC 1108 MATH-PACK subroutine JACMX [ 3] . This subroutine com- 
putes the eigenvalues of the matrix of inertial properties; further, the corre- 
sponding eigenvectors are precisely the direction cosines (between the orbiting 
cqordinates of Figure 3 and the principal axes) required in the gravity gradient 
torque computations to follow. It is noted that in the inertial hold mode, the 
inertia tensor must be transformed, to be in accordance with this procedure, 
to the orbiting coordinates; however, the Y-POP mode’ s tensor is already 
computed with respect to the orbiting T reference frame. 

The method used by JACMX is a modified Jacobian method in which the 
matrix is diagonalized. This is accomplished by using elementary orthogonal 
transformations to annihilate successive off-diagonal elements. ,The matrix 
of eigenvectors is the product of all the transformation matrices. 

GRA V ITY G RA D I ENT TORQ UE 


It is well known that the gravity gradient torque acting on a rigid bocfy 
is solely a function of its mass distribution, once its position in space is known. 
One commonly used form for this torque [ 4] is 



in which, in the nomenclature of the present work, 


W = orbital speed of the space station for a circular orbit, 
which is 27 t over the orbital period t ; 


(i. J» k) = 


unit vectors in the T-frame of Figure 3, which are the 
same as those of frame V for the Y-POP mode; 


(I, J, ^ 


unit vectors in the principal directions; 


(Ii» h, I3) 


principal moments of inertia with respect to the mass 
center, about (I, J, K), respectively; 


M 


c 


gravity gradient torque vector with respect to the 
body’ s mass center. 


Nothing more need be said here except that the computer is programed 
to automatically compute the eight required dot products in equation (37) as 
soon as the eigenvalues (principal moments of inertia) and eigenvectors 
(principal directions) are found via the library program described in the pre- 
ceding section. 

After performing the arithmetic required in equation (37), the computer 
program yields the two gravity gradient torque components and their resultant. 

Following the printing of these results, the program then increments 
d and repeats the procedure until an orbit of Y-POP or inertial hold has been 
completed. The results obtained for several examples will now be presented 
and discussed. 


RESULTS, CONCLUSIONS, AND RECOMMENDATIONS 

The first computer run made was for a space station consisting only of 
a power module and solar panels. The angles were set simply to ^2 ” ~ ^y 

= \ = 0 , thereby placing ever 5 dhing in the same plane for convenient checking. 



The program was then debugged successfully by checking the masses, inertial 
properties, and gimbal angles over Y-POP and inertial hold orbits. Then the 
parameters were changed to (f>^ = 23. 5 deg , A = 35 deg , T _ = 90 deg , and 


A ^ = 35 deg 


r = 90 deg 

y ® 

Ay = 0 , and the results for the same station indicated a rise in gravity torque 


magnitude from zero at 0 = 0 to as much as 0. 19 N-m in Y-POP and 0. 28 N-m 
in inertial hold. The Y-POP change is caused solely by solar panel gimbaling, 
showing its importance for symmetric stations. 


The space stations sketched in Figure 6 were then run for two reasons: 

(1) to further insure that the procedure and the entire program work correctly, 
and (2) to investigate such effects as solar panel gimbaling and asymmetry of 
the configuration in station buildups of practical importance. 

The magnitudes of the gravity gradient torque on No. 12 for Y-POP and 
inertial hold modes are shown in Figure 7, along with the orbital parameters 
used and other data of interest. It is immediately observed again from the 
Y-POP curve that solar panel gimbaling has a definite effect on the torque. 

This follows from the observation that the station in the Y-POP mode would feel 
the gravity torque value at 0 = 0 throughout the orbit if the panels were not 
allowed to gimbal. However, with gimbaling, there is a spread in gravity torque 
magnitude of about 10 percent over the entire orbit of station No. 12. It is 
observed that the gimbaling has an even greater effect on station No. 18 
(Fig. 8) , the spread in torque being about ±15 percent over its orbit. Note 
from Figure 9 that the largest changes in Y-POP torque correspond closely to 
the largest changes in the gimbal angles, as expected. Hence, we conclude that 
solar panel gimbaling is a significant factor to be considered in sizing control 
nioment gyros to meet momentvun requirements. 

Figures 7 and 8 also indicate interesting results for the inertial hold 
modes for these two stations. (Note that both curves indicate the periodicity 
with respect to tt that was expected. The torque should indeed be periodic, 
with twice orbital frequency. ) A very important result that can be seen in 
these figures is that there is about a threefold increase in the bias, or noncyclic, 
portion of the torque when going from station No. 12 to station No. 18. A run 
made for an even more unsymmetrical station (No. 47, not shown) showed an 
even higher bias torque than No. 18, having peak-to-peak magnitudes of 0. 15 
to 11. 52 N-m. A run was also made in which stations Nos. 12 and 18 were 
locked into inertial hold modes at another orbital angle, with similar results. 
Thus, a reasonable conclusion seems to be that each proposed asymmetric sta- 
tion configuration should be separately examined and that "average geometry" 
cases should be avoided. 
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Figure 


6. ISS configurations 


studied. 







Figure 8. Gravity gradient torques on No. 18. 

Some recommendations for extensions of this work are now tabulated. 
In future work, we wish to include: 


1. The capability to trim the axes 

2. The computation of angular momentum requirements 

3. The ability to feather the solar panels so as to minimize angular 
momentum buildup during loss of sunlight 
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, 081 - 


4 . The forming of the GSS station by adding a second power module 
and solar panels 

5 . The capability to directly feed ^the beta angle (77 ) instead of the 

parameters used by the writer; also to feed ^ a precalculated inertia tensor 
or to feed out the internally computed inertial properties and/or gimbal angles 
for other uses 

6. The writing of the differential equations of rigid body motion of the 
station (possibly including aerodynamic torques) and their examination for 
stability. 

Most importantly, we wish to be able to rotate the station through any 
angle at any point in its orbit and have it remain fixed in that position for a 
given time, during which torque and momentum would be computed. 
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